CERN-TH/99-414 
hep-ph/000115^ 
December 1999 



(N 



O 
O 



Applications of QCD|^ 



o 

o , 

O ! Martin Beneke 



Theory Division, CERN 
CH-1211 Geneva 23, Switzerland 



^ . Introduction 

■ Quantum chromodynamics (QCD), the theory of the strong interaction, attracts a 

large body of theoretical work and experimental investigation. "Applications of QCD" 
O ■ [theory] probably means that the former should have something to do with the latter. 

This is a severe restriction. It leaves out, for example, phenomena at high temperature 
and large baryon density, which have not been tested in terrestrial laboratories. Yet 



p i' there have been fascinating developments in these directions, which reflect once more 

p the complexity that follows from a Lagrangian as simple as —(7^/4 + quarks. In 

D . this talk I concentrate on high energy QCD processes, which means that at least 

^ ' some part of the process must be tractable perturbatively. Even within this narrow 

frame striving at completeness would do injustice to the diversity of the (sub-)field. 
^ ■ The following gives a survey and assessment of recent theoretical results on selected 

■ topics. For details please consult original references and topical reviews. Apologies 

for omitting topics that should have been included, but have not been for various 
reasons (lack of time, competence, ...). 

The conceptual basis for discussing QCD processes at large momentum transfer 
Q is provided by factorization: 

da = da{Q, jj) ® F{fi, Aqcd) + 0{Aqcd/Q). (1) 

The first factor, dd-{Q,fi), is insensitive to long distances of order 1/Aqcd- It is com- 
puted in perturbation theory as scattering of quarks and gluons and depends only on 
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the strong coupling ag and heavy quark masses. The second factor accounts for the 
fact that experiments are prepared and measurements are done far away 1 fm) 
from the interaction point. Aqcd) parametrizes this long-distance sensitivity 

in terms of process dependent quantities: vacuum condensates, parton distributions, 
fragmentation functions, light cone wave-functions and many more. In principle, 
F(yU, Aqcd) depends only on and light quark masses, but since ^^(Aqcd) is large, 
we cannot compute it in perturbation theory. However, being independent of the hard 
scattering process, the same Aqcd) may appear in a generic class of processes. 
There are some fortunate cases in which long distance sensitivity appears only in 
0{A/Q) in (H). In these cases, we have particularly clean predictions, if Q is large 
enough. In general, we need to provide F{fi, Aqcd). It can sometimes be computed 
non-perturbatively by numerical methods ("lattice QCD"). Or it may be approxi- 
mated by models of low-energy QCD. More often, however, some measurements are 
used to determine Aqcd); others are then predicted. This makes QCD seem 
to depend on many infrared parameters along with a^. It also implies iterations of 
theory and experiment to arrive at predictions. 

Eq. (0) suggests a procedure: for any given large momentum transfer process (i) es- 
tablish (|l]), identify Aqcd); (h) compute da{Q,fi) accurately; (iii) if Aqcd) 
is known, predict da, otherwise determine Aqcd), if da is measured; (iv) check 
the accuracy of this procedure by addressing power corrections 0(Aqcd/<5)- The 
outline of this talk is divided in sections according to this procedure rather than by 
topics, although in different order. Sect. [l| covers perturbative calculations. Sect. |^ 
power corrections. Perturbative expansions of da{Q, /i) often fail in special kinematic 
regions, but accurate results can be recovered upon all-order resummations. In Sect. |^ 
I discuss three representative examples of this situation. Finally, Sect. ^ is devoted 
to some processes for which factorization has been established more recently. 

It is important to remind ourselves that working with QCD we take many things 
for granted which have never been proven: that (|1]), obtained from factorization 
properties of Feynman diagrams, holds non-perturbatively; that the operator product 
expansion holds non-perturbatively; that perturbative expansions are asymptotic; 
that lattice QCD approaches the correct continuum limit. The overall picture of 
consistency that has emerged in applications of QCD suggests a pragmatic attitude 
towards these problems. However, the questions remain. 

1 Perturbative calculations 

For long-established QCD processes there are no easy perturbative calculations any 
more. Increasing the accuracy by one order in as has become technically demand- 
ing, usually requiring extensive or automated algebraic manipulations by computers 
and/or numerical computing. The complications increase by increasing the number 
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of loops, or the number of mass scales or external legs. 



1.1 More loops 

Totally inclusive quantities are related to imaginary parts of correlation functions. 
This avoids infrared divergences in intermediate expressions. Such quantities are 
candidates for fully automated evaluation 

The correction to e^e~ hadrons (massless quarks and gluons) and related 
observables, and to some deep inelastic scattering sum rules, have been known for 
some time 0. More recently, the QCD /5-function P] and quark mass anomalous 
dimension |Q have been computed at 4-loop order. 

These results use that any 3-loop, massless, 2-point integral is calculable in dimen- 
sional regularization. The most important tools are the integration-by-parts method 
[^], infrared re- arrangement which reduces the calculation of the 4-loop pole part 
to the above class of diagrams, and powerful computers that handle the algebra con- 
nected with about 10^ Feynman diagrams. Another important class of diagrams 
which is generically calculable is 3-loop, massive, vacuum bubble diagrams 0. There 
is no obvious way to extend these results to one more loop. 



1.2 More scales 

Observables that depend on more than one kinematic invariant or a kinematic invari- 
ant and quark masses are difficult, even if they are totally inclusive. A method that 
has led to a number of interesting new results is based on asymptotic expansions in a 
ratio of scales, such that each term in the expansion is a single-scale integral that is 
analytically solvable. This method can be used even if the expansion parameter is not 
small if many terms in the expansion can be obtained and if the radius of convergence 
is sufficiently large or convergence can be improved by Fade approximants. 

Asymptotic expansions can be performed (i) for large external momenta, small 
masses or for large masses, small external momenta M; (ii) around mass shell P; (iii) 



near thresholds |10] or in t/s for 2-^2 scattering; (iv) in Sudakov limits |11|. These 
expansions are done on the integrand level. The fact that loop momenta cover all 
scales implies that, in general, extra terms have to be added to the Taylor expansion 
of the integrand. 

A nice example to illustrate the method is the 3-loop coefficient in the relation 
between the pole mass and the MS mass of a heavy quark [|l^. This requires 3- 
loop on-shell integrals, which are not known. Instead expand the quark self-energy 
around external momentum = 0, which reduces the problem to 3-loop vacuum 
bubbles, which are calculable. Then put p'^/m^ = 1 and use Fade approximants. 
Expansion to order {p^/m^Y^ (plus information from the opposite limit 3> m^) 
gives r3 = 3.10±0.06 for the coefficient at order a1 {rif = 4). (In retrospect, this turns 
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out to be a "bad" example, because the 3-loop on-shell integrals are in fact exactly 
calculable [|l^, |14|. The exact number is = 3.0451 . . . [|l^], in nice agreement with 
the previous semi-analytic result.) 

Applications of this method up to now concern quantities with internal masses or 
on-shell lines. e^e~ bbX has been obtained at order for general and in 

an expansion near threshold . The corrections to inclusive heavy quark decays 

111, t ^ 6iy PI and 6 ^ m/z/ @. The last result is 



clu 



has been calculated for b 
particularly impressive, because the asymptotic expansion is obtained algebraically 
to all orders. It is then resummed to an exact result. 



1.3 More legs 



Higher order jet calculations pose a different sort of challenge, because the kinematics 
becomes complicated (as the number of jets increases), and because the calculation is 
done on the amplitude level. Infrared singularities cancel in an intricate way, or are 
factorized into parton densities (fragmentation functions) after cancellations. Almost 
certainly the final result is obtained after numerical integration. 

Relatively recent results include NLO corrections to e+e~ — > 4 jets 



e ' e 



jets with quark mass effects [^, which provide us with a first, yet imprecise, evidence 
of scale-dependence of the bottom quark mass Partial NLO results exist on 

pp ^ 3 jets The full result is supposed to be completed soon. 



1.4 Towards NNLO jets 

The conceptual and technical frontier is set by NNLO jet calculations, the basic 
process being 2 — > 2 {pp —>■ 2 jets or 1 jet inclusive, pp —>■ 77X) or 1 ^ 3 (e^e^ 3 
jets). NNLO calculations provide detailed insight into jet structure and a better 
determination of a^- In e~^e~ —>■ 3 jets they are important to understand the interplay 
between perturbative and power corrections. 

There are several components to the NNLO jet project. The amplitudes have to 
be computed, which include 2-loop 4-point diagrams. Amplitudes with five and six 
partons have to be integrated analytically over the singular regions of phase space. 
After cancellation of infrared divergences, the remaining phase space integrals have 
to be evaluated numerically efficiently. There has been progress on many of these 
components recently. 

Because of the integration over singular regions of phase space even tree ampli- 
tudes are non-trivial. In 2 — > 4 tree amplitudes one encounters a new situation, when 
two partons become simultaneously soft or three partons become coUinear. The last 
case (squared and integrated over phase space) gives rise to a new class of sphtting 
amplitudes (functions), when one parton decays into three coUinear partons, which 
generalize the usual splitting functions. All soft, coUinear, and mixed, limits have now 
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been analyzed |2^. Likewise, although the 1-loop five-point amplitudes are known 
in four dimensions, this is insufficient, because the two jet cross sections includes 
configurations where two of the partons are not resolved. Making use of the univer- 
sality of soft and coUinear limits, these amplitudes are now known to all orders in the 
dimensional regularization parameter e in those kinematic regions, where the phase 
space integration is singular . 



The most difficult amplitude is the 2-loop virtual correction to the basic 2^2 
(or 1^3) process. Until very recently, it has been unclear whether the basic scalar 
double box integrals are analytically calculable. In a stunning calculation an 



analytic result was obtained for the planar double box integral, expressed in terms of 
elementary special functions, and an algorithm was provided to compute the integral 
with arbitrary numerator |^^. It is equally surprising that this result was obtained 
by elementary methods: the a-representation and Mellin-Barnes transformation and 
summations of multiple sums obtained after taking the Mellin-Barnes integrations. 
The crossed double box was subsequently calculated using the same methods. 



The numerator algebra, connected to multiple products of three-gluon vertices, is, 
however, highly non-trivial, and remains to be done. Methods, based on helicity 
amplitudes, colour decomposition, special gauges and unitarity exist to simplify the 
task [^. Up to now this has been completed in a toy = 4 supersymmetric theory 



(leaving the scalar integral unevaluated) [^, and more recently for the maximally 
helicity violating amplitude in QCD pT[| . 

Many of these results can be used also for NNLO corrections to e~^e~ 3 jets. 
However, the 2-loop double box integrals with one off-shell external leg are not yet 
known. The infrared singularities at order are known ||3^, but the structure 

of 1/e poles remains to be elucidated. 

In my opinion, the results that have been achieved over the past two years make 
success predictable, at least for NNLO 2 jets. On the other hand, many hard algebraic 
and numerical tasks remain to be done. Even with a concerted effort the relevant time 
scale is years rather than months. However, this is clearly a beautiful case, where 
most of the most advanced techniques for perturbative QCD calculations merge into 
a single project. 

1.5 NNLO parton evolution 

NNLO jets require NNLO parton distributions. Evolution of these parton distribu- 
tions requires the NNLO DGLAP splitting functions. The complete NNLO splitting 
functions are still unknown, although some moments have been computed p3| some 
time ago and further constraints exist in the large-a; and small-x limit. 

Large evolution means large Q^, since parton distributions are typically deter- 
mined experimentally at some low scale. Large means large x and hence the 
known moments may already provide accurate information. Indeed, first construc- 
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tions of approximate NNLO non-singlet splitting functions have appeared that make 
use of the known information |3^. The constraints at small x turn out to be quite 
weak, but there seems to be little uncertainty for x > 0.1. When the splitting function 
is folded with a typical parton distribution this range increases to x > 0.02. 



2 Beyond leading power 

Perturbative expansions, if computed to arbitrarily high order, ultimately diverge. 
They become useless beyond a certain order, unless they are summed. In recent 
years, we have learnt to turn this embarrassment into a benefit, since the pattern 
of divergence tells us something about the scaling of power corrections (Aqcd/Q)" 
to a hard scattering cross section. A particularly interesting type of divergence, 
called infrared renormalon, is related to integration over small loop momentum in 
Feynman integrals |3^. Roughly speaking, there is a relation between perturbative 



long distance sensitivity, the size of perturbative coefficients in higher orders, and 



the scaling of non-perturbative power corrections For inclusive deep inelastic 

scattering, n = 2, and one recovers higher twist corrections predicted by the operator 
product expansion. 

2.1 Event shape observables and energy flow 

For other, less inclusive, observables, such as event shape variables in e^e~ and ep 



collisions, one often finds n = 1 |3^, Since these variables are order as per- 



turbatively, they are prone to large non-perturbative (and perturbative) corrections. 
They have been investigated intensively over the past two years, theoretically and 
experimentally. 

The leading power correction originates from soft partons emitted from a fast, 
nearly back-to-back qq pair. Write 



(5) = {S^-Hf^i)) + I {S^'^'if^i)) + 0(A^cd/Q') (2) 



for an average event shape variable S. The experimentally measured energy depen- 
dence of (S) clearly supports the existence of a 1/Q power correction with a reasonably 
sized normalization {S'^^{fij)), which is non-perturbative. An interesting hypothesis 
(also applied to event shape distributions states that the non-perturbative cor- 



rections are universal, i.e. {S'^^{fj,j)) oc csa{fJ.i), where cs is observable dependent, 
but calculable and a{ni) is non-perturbative but independent of S . (This applies 
to thrust, jet masses and the C parameter. Other event shapes, such as jet broaden- 
ings, involve complications [0].) Four parton final states with two soft partons have 
also been investigated [|^, ^ Remarkably, one finds that cs is rescaled by the 



same factor for a variety of shape observables P3, WM. The universality hypothesis 
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has led to a number of instructive experimental tests. Recent results on average event 



shapes and event shape distributions in e~^e~ annihilation |^ and DIS [^] tend to 
confirm the hypothesis within the expected accuracy. However, the fact that the 
value of as, fitted simultaneously to each S, is somewhat unstable indicates that the 
present understanding is not perfect. 

Universality may hold for a special class of observables, but it would be surprising, 
if it held in general. What is needed to shed light on the issue is a factorization 
theorem for soft gluons beyond leading power. Recall that factorization theorems 
for event shapes usually demonstrate that soft gluon corrections cancel at leading 
power. We are now interested in the leading contribution that is left over after this 
cancellation. 



In the problem is approached in terms of energy flow of soft particles. The 



universal, non-perturbative objects relevant to the two-jet limit {qq plus soft partons) 
are 

k 

G{n,, ...,nk; fn) = {0\W^ JJ ^(^«) ^|0)' (3) 

i=l 

where measures soft energy flow in the direction of n^, fij is a factorization 

scale that defines what "soft" means, and W denotes a product of eikonal lines for 
the energetic qq pair. The G{ni, . . . ,nk] fii) are horribly complicated objects and it is 
hardly conceivable that they could ever be extracted from measurements. However, 
the fact that they are independent of the hard scale Q already entails interesting 
predictions. For example, event shape distributions can be expressed as a convolution 
of a perturbative distribution and a non-perturbative Q- independent, but observable- 
dependent "shape function" , that follows from these energy flow correlation functions. 
Event shape averages can be represented as 



(5)i/Q = J dnws{n)G{n), (4) 

with a calculable weight function ws{n). The single energy flow correlation function 
G{n) can in principle be determined from the leading power correction to the energy- 
energy correlation. 

I find this a promising step towards understanding soft power corrections. The 
concept of energy flow is clearly important and deserves more attention, as it corre- 
sponds directly to calorimetric measurements. Observables that can be represented 
in terms of energy flow are automatically infrared safe. They may also be defined 
non-perturbatively and therefore be amenable to a more systematic analysis of power 
corrections 



2.2 OPE of the plaquette 

There are also things that don't work as expected. Consider the operator product 
expansion (OPE) of the plaquette expectation value in pure gauge theory at finite 
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lattice spacing, i.e. the inverse lattice spacing takes the role of the scale Q ^ Aqcd- 
The OPE gives 



(plaquette) = ^ c^af ^Q)" + {-GG) + .... (5) 

n=l Q ^ 

The coefficients c„ of the perturbative expansion have been computed to 8th or- 
der numerically [^]. After transformation to a continuum coupling definition, the 



coefficients exhibit the expected infrared renormalon growth. Summing the series 
approximately should give an accuracy of order hence subtracting the summed 

series from non-perturbative Monte Carlo data for (plaquette) , the remainder should 
scale as l/Q^, consistent with the scaling of the gluon condensate term. 

Contrary to this expectation, the remainder is found to approach a perfect l/Q"^ 
scaling behaviour . Since the OPE is one of the few tools we have to go beyond 



perturbation theory, this is clearly something we should understand. There may be 
subtleties with the transformation to the continuum scheme, since this transformation 
is not known to 8th order or may also have power corrections. The effective action 
at finite lattice spacing contains an infinite set of higher dimension operators. Could 
these add up to a l/Q^ power correction ||5T[ so that the result is a lattice artefact? 



But there may be less profane explanations such as power corrections from short dis- 
tances that affect coefficient functions (and therefore would not contradict the OPE) 



3^ . This possibility is not ruled out by any argument. It presents a fundamental 
question that challenges our understanding of non-perturbative short- distance expan- 
sions. It would also have implications for the phenomenology of power corrections to 
current correlation functions. For these reasons, the problem raised by |5D[ should be 
cleared up! 



3 Perturbative resummations 

Returning to perturbative expansions, it is not unusual that a perturbative expan- 
sion in as breaks down, even though the coupling constant is small. This happens 
because the smallness of the coupling constant is compensated by a large kinematic 
invariant. In effect, one is dealing with a multi-scale problem. If all scales are large 
compared to Aqcd, the problem is perturbative and may be subjected to system- 
atic all-order resummations. The kinematic conditions leading to the breakdown of 
perturbation theory can be quite different and the resummations reflect completely 
different physics. In this section I discuss three examples of such resummations, where 
progress has been made over the past two years. 
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3.1 Parton thresholds 



A familiar source of large kinematic corrections is related to partonic thresholds. 
Consider the differential cross section 

= J2 fi/^ ® fj/B ® d&ij^f (6) 

for a hard hadron-hadron collision. Large logarithms appear in da, when the cms 
energy s of z + j is just large enough to produce a given final state. For exam- 
ple, in production of a massive vector boson with mass Q, the leading correction is 
ag(ln^"'~^(l — -2))/[l — z]^ at order a", where z = Q'^/s, and perturbation theory 
breaks down for z ^ 1. 

In this case large logarithms originate from the lack of phase space for real emis- 
sion and the incomplete cancellation of sensitivity to collinear and soft momentum. 
Because of this relation the structure of these logarithms is well understood. The 
logarithms exponentiate and can be resummed: 

Jdzz^-Ua{z) = 

H{as) exp[\nN g^iaJnN) + g^iaJjiN) + aMc^s'^nN) + ...]+ 0(1/ N). (7) 

This resummation was worked out at next-to-leading logarithmic order (i.e. including 
g2{c(slT^ N)) some years ago for 2 —>■ 1 processes (massive vector boson production) 
and 1 —>■ 2 processes (event shape variables in e~^e~ in the 2-jet limit) ||54|| . 
Next-to-leading logarithmic resummation has now been extended to 2 — > 2 scat- 
tering processes [0, |56[. Several new complications appear in this case. Since the 



underlying hard process depends on an additional kinematic invariant, (— t)/s, so do 
the functions that appear in the exponent of (0). Furthermore, the 2^2 ampli- 
tude contains several colour amplitudes and since soft gluon emission carries away 
colour, these amplitudes mix, turning the exponential into a matrix exponential on 
the independent colour amplitudes. While the structure of resummation remains thus 
the same, the technical complications make the formalism more difficult to apply in 
practice. 

Fortunately, simplifications occur for total cross sections. NLL resummed results 
have been presented for heavy quark production [^, ^| and prompt photon produc- 



tion |5^. For di-jet production at large transverse momentum, the formalism is in 
principle complete, but it has not yet been implemented It turns out that at 

energies of interest for heavy quark production and prompt photons, the effect of re- 
summation is typically small, i.e. within the renormalization scale variation of a fixed 
order NLO calculation. The real benefit of resummation is a significant reduction of 
this scale dependence compared to NLO QCD, and hence, probably, the theoretical 
uncertainty. The Ej- spectrum of prompt photons at low remains in disagreement 



9 



with the data |6^. Since at Et ~ several GeV power corrections in 1/Et, or intrinsic 
transverse momentum, can be very important, this is hardly a serious issue. It is, 
however, a serious problem for determining the gluon distribution at large x. 



3.2 Non-relativistic 

A different kind of partonic threshold is encountered in heavy quark production in 
e~^e~ annihilation. When the cms energy is just larger than 4mg, the quark and 
antiquark move at small relative velocity and attract each other through a strong 
Coulomb force, even if as is small. Formulated as a perturbative resummation prob- 
lem, we need the terms 

Re+e-^QQX ~ ^ E " " {l(LO); t;(NLO); al a^v, v^NNLO); . . .} (8) 

n=0 \ ^ 

at leading order (LO), next-to- leading order (NLO), etc., where v is the small relative 
velocity. 

The LO resummation is done by solving for the Green function of the Schrodinger 
equation with the Coulomb potential. To be more systematic, such concepts from 
quantum mechanics have to be derived from QCD, incorporating correctly the short- 
distance structure of QCD. This is done by a sequence of non-relativistic effective 
field theories. Quarks and gluons can be classified as hard (h), soft (s), potential 
(p) and ultrasoft (us) |T0|. Then these modes are integrated out successively, ac- 
cording to the scheme CQc-oiQih, s,p); g{h, s,p,us)] £nrqcd[Q(s,p); ms)] 
'^p'NRQCDlQip)', g{us)], passing from QCD to non-relativistic QCD to potential- 
non-relativistic QCD The equation of motion of PNRQCD is exactly the Schro- 
dinger equation, with corrections to it that encode the information of the short- 
distance modes that have been integrated out. 

With the help of this method the NNLO resummation has been performed. This 
leads to first principle NNLO calculations of ti production near threshold (in e+e~ 
collisions) |6^, |6^. The NNLO correction turns out to be very important and 

has led to the conclusion that it is the MS top quark mass rather than the pole mass 
than can be determined more accurately, though indirectly, from the cross section 



near threshold p4|, |65[ . Another important application concerns the determination of 



the b quark mass from e~^e bbX [67, The recent analyses that take care 



of adequate bottom mass renormalization prescriptions converge towards a common 
value for the bottom quark MS mass, which I average as mb{fnb) = 4.23 ± 0.08 GeV. 
The centre of attention is now on understanding logarithmic corrections in v [^, |69| . 
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3.3 High energy, small x 

The high energy hmit oi QCD cross sections is an old, yet unsolved problem. 

Large logarithms can appear either in the high-energy limit of hard partonic reactions, 
such as in 7*7* scattering or forward jet production, or in the small-x behaviour 
of parton distributions and their evolution. The leading logarithms (aglns/Q^)" 
have been summed long ago by means of the BFKL equation This leads to 

cross sections that rise as s"""^^^"^ [dg = Ncag/Tr) with energy. For many years most 
theoretical work has been concerned with the physical mechanism that would make 
the high energy limit compatible with unitarity, but a quantitative theory has not yet 
emerged. Most of the recent activity in small-x physics has however been inspired by 



the completion of the NLO correction to the BFKL kernel |71|, and its interpretation. 
The following discussion concentrates on this aspect. 

Recall that phenomenological applications of LO BFKL theory have remained 
ambiguous or unsuccessful. HERA data on the gluon density indicates that DGLAP 
evolution works well, in fact too well, down to x ~ 10~^. No resummation of Inx 
corrections to the evolution kernels is required. There is some flexibility in the input 
gluon distribution, nevertheless the message is that departures from DGLAP cannot 
be large. Virtual photon scattering has been measured at LEP [^. Even allowing for 



the fact that LO BFKL may not predict the normalization of the cross section well, 
the observed energy dependence is less steep than predicted. Forward pion production 
at HERA [|73| may be described by LO BFKL, but other interpretations of the data 
seem possible. 

It is therefore clearly interesting to see how NLO corrections affect this compari- 
son. In the high energy limit, the cross section factorizes schematically as 



where A usually represents a virtual photon and B a virtual photon or a proton. 
In the latter case the impact factor $p(/c2) is not perturbatively calculable. ki^2 
denote transverse momenta of the scattering objects, ~ Q for virtual photons, 
k ~ Aqcd for protons. The factorized form (|^) is believed to hold to next-to-leading 
logarithmic order, but beyond this order there are terms that cannot be associated 
with the four-(reggeized)-gluon Green function G^{ki,k2). G^{ki,k2) satisfies the 
BFKL equation 

ioG^ku k2) = 6^^\k, -k2)+ f —K^k^, k) GUK k2). (10) 

J 71 

Roughly speaking, the leading order kernel K^{ki,k) sums a single gluon ladder 
exchanged between A and B with emissions ordered in longitudinal momentum. The 
NLO correction has to account for all configurations in which one power of Inx is lost. 
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Partial results have been collected over many years and the full NLO correction has 
finally been completed [^]. It is usually presented through the action of the kernel 



on a set of test functions: 

J d'k'KUk,k') (l^j =asXo{l) 



^2' 

1 - PoOis In — 



+ aixi{l). (11) 



In the saddle point approximation for the inverse Mellin integrals, treating as 
small, the energy growth of hard high energy cross sections is then determined by 

A = a,xo(l/2) + «,V(l/2) = 4 In 2 [1 - 6.5«,] , (12) 

ignoring the scale dependent part of the kernel. The NLO correction is huge, large 
enough to modify qualitatively the conclusions drawn from leading order, which is 
good. At the same time, the NLO kernel taken at face value leads to non-sense results 



[[7^ , unless cxg < 0.05, which is unrealistically small. 

Much effort has gone into the question whether the NLO result invalidates the 
BFKL resummation programme as a whole. To answer this question one has to go 
beyond a systematic resummation of high energy logarithms. Such a step is unavoid- 
ably ambiguous and needs to be motivated by physics arguments. It appears that 
much of the NLO characteristic function Xi(7); even near 7 = 1/2, can be under- 
stood from the singularities at 7 = and 1. Note from (^) that these singularities 
correspond to transverse logarithms. The leading singularities 1/7^'^, 1/(1 —7)^'^ are 
related to the symmetric energy scale ^1^2 chosen in (^, the running coupling and 
the non-singular terms of the LO DGLAP splitting function. Remarkably, Xiil) is 
extremely well reproduced just by keeping these singularities. 

This suggests that these singularities should be summed to all orders. Unphysical 
transverse logarithms generated by the symmetric energy scale can be removed 
by replacing 

Xoil) ^ Xo (7) = 2V^(1) - V^(7 + uj/2) - V;(l - 7 + uj/2). (13) 

Although not unique, this seems to be a particularly natural choice. After performing 
this replacement, the NLO correction is reduced, though not small. Further support 
for this resummation arises from the possibility to introduce a "rapidity veto" — 



Ui > A which is essentially a hard cut-off on the momentum region, where 
the ordering in rapidity was not a good approximation in the first place. After 
resummation, the A-dependence is small and the NLO correction moderate for all A 



7^ , which indicates that the resummed kernel is less sensitive to momentum regions 
where the approximations necessary to derive it are not valid. 

The remaining 7-singularities are double poles. Two further modifications beyond 
NLO small- a; logarithms need to be performed to take care of them. First, rather than 
improving the DGLAP anomalous dimension by small-x logarithms, we can take 
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the opposite point of view and improve xil) by taking into account all information 
on collinear logarithms ||78|, |7^, |8^. In this way one can arrange, in addition, for 



momentum conservation, which requires vanishing anomalous dimension at a; = 1. 
Second, the 1-loop evolution of ag can be taken into account exactly, rather than 
perturbatively as in (0). There are two cases to consider, symmetric processes with 



fci ~ ^2 ^ Aqcd eM asymmetric processes |]82|, |79[. In the latter case, with 



Q ~ ^1 ^ ^2 ~ Aqcd as for deep inelastic scattering, one must also apply collinear 



factorization to the four gluon Green function, such that [33 

G^ih, h) = FUV(fci) ■ (fc2) + Oikl/kl). (14) 

The dependence on the non-perturbative low momentum evolution of the running 
coupling is factorized into F^{k2), which can be absorbed into the input gluon dis- 
tribution. This part remains beyond perturbative control, although it may well con- 
trol the actual small-x behaviour of the gluon distribution. On the other hand, only 
F^^{ki) is Q-dependent and hence determines the evolution of the gluon density. 
There seems yet not to be an unanimous opinion on which of these aspects is 



most important. For example, argues that A should be considered as a non- 
perturbative parameter, while [|79| takes a less agnostic attitude. Ref. |8^, on the 
other hand, emphasizes the role of the running coupling, demonstrating that the 
effective scale of in the anomalous dimension increases as x decreases, because of 
ultraviolet diffusion. It is also claimed that this leads to an improved fit to structure 
function data compared to a standard DGLAP fit, which is definitely interesting. 
Despite these different viewpoints, theory clearly seems to be on the right track, as 
the results consistently point towards a smaller (but positive) hard pomeron intercept 
compared to LO BFKL. The resummed gluon anomalous dimension is also close to the 
DGLAP one down to rather small moments. It will be interesting to see consolidation 
of this field and the first true NLO+improved BFKL predictions for physical processes 
(which needs as yet unknown NLO impact factors). 



4 Novel factorization "theorems" 

In the past sections I discussed hard scattering processes which have been known as 
such. But for other processes factorization of its short-distance part has been estab- 
lished only recently. Often factorization comes at the expense of introducing new 
non-perturbative parameters. Even if these parameters are not accessible immedi- 
ately, much is gained in terms of conceptual clarity. In this section I discuss three 
examples of such "new" applications of QCD. 
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4.1 Hard diffraction 



A particularly nice example is hard diffraction ||8^. Discovered in hadron-hadron 



collisions by UA8 about a decade ago ||8^, after the inspiring work of the extent 



to which hard diffraction is a hard process, has remained rather unclear. This has 
changed completely with the arrival of accurate data on hard diffraction in ep scat- 
tering the demise of Regge terminology, and the realization that hard diffraction 
in DIS can be described in close analogy with inclusive DIS [^, 

In hard diffractive DIS, 7*p — » Xp, the proton scatters (quasi-)elastically off 
a virtual photon, which fragments into a colour neutral cluster X. The scattered 
proton is usually not detected, but since it typically loses only a small fraction of its 
momentum, the event is identified by a large gap in rapidity between p and X. About 
10% of all DIS events are rapidity gap events. Furthermore, hard diffraction is not 
suppressed with l/Q"^ relative to inclusive DIS. 

In close analogy with inclusive DIS, the diffractive cross section factorizes into a 
short- distance cross section and a diffractive parton distribution |^9|, pO| : 



E/%*-(Q.x.,;,)*iM^. (15) 

The diffractive parton distribution //'(l/, ^, /i) represents the probability to find 
parton i in the proton with momentum fraction y under the condition that the proton 
stays intact and loses longitudinal momentum fraction ^. Note that this definition 
makes no reference to Regge factorization or the pomeron. Neither does it make 
reference to a rapidity gap Ay , which follows from kinematics alone when ^ is small: 
Ay ~ ln(l/^). The hard scattering occurs on a single parton as in ordinary DIS. 
The dynamics that is responsible for the formation of a colour-singlet cluster is non- 
perturbative and therefore part of the definition of the diffractive parton distribution. 

The physical picture of hard diffraction is perhaps clearest in the proton rest frame 
and reminiscent of the "aligned jet model" pT] , p2| . In the proton rest frame, at small 
Bjorken x, the virtual photon splits into a qq pair long before it hits the proton. The 
qq wave-function of the virtual photon suppresses configurations in which one of the 
quarks carries almost all momentum. Yet it is these configurations that give rise to a 
large diffractive cross section, because the wave-function suppression is compensated 
by the large cross section for the scattering of a qq pair of hadronic transverse size 
off the proton. The harder of the two quarks is essentially a spectator to diffractive 
scattering. The scattering of the softer quark off the proton is non-perturbative 
and cannot be described by exchange of a finite number of gluons. Hence there is an 
unsuppressed probability that the softer quark leaves the proton intact. This explains 
the leading twist nature of hard diffraction. The details of the scattering of the softer 
quark off the proton are encoded in the diffractive quark distribution. In a similar 
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way, the qqg configuration in the virtual photon, in which the qq pair carries almost 
all momentum, gives rise to the diffractive gluon distribution. 

Because the short-distance cross section a"^ * of hard diffractive DIS is identical to 
inclusive DIS, the evolution of the diffractive parton distributions is identical to those 
of ordinary parton distributions. It follows that the characteristics of diffraction 
are entirely contained in the input distributions at a given scale. It is therefore 
interesting to model these distributions. The original idea of a partonic content 
of the pomeron 
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can be interpreted as an ansatz in which the diffractive parton 
distribution factorizes into a pomeron flux factor, which determines the ^ dependence, 
and a parton distribution in the pomeron which depends only on f3 = x/^. The precise 
data from HERA do not support this simple ansatz any more, although the problem 
can be fixed by adding more Regge poles. More recent approaches model the proton 
field off which the Fock states of the virtual photons scatter. The semi-classical 
approach which preceded the factorization theorem, can be formulated in such 
a way that it models the diffractive parton distributions It can be justified 

for a large nucleus [03. Applied to the proton it gives a reasonable description of 
both diffractive and inclusive DIS [Q. (See [0 for earlier work that contains some 
elements of the semi-classical approach.) Another approach is based on two gluon 
In this case one either has to deal with an infrared divergence. 

Remarkably, these three 



exchange 
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or couple the gluons to a small size toy nucleon as in 
approaches give similar results on the /3-dependence of diffractive parton distributions 
and agree on the fact that the gluon distribution is enhanced by a large colour factor. 
This leads to positive scaling violations already at relatively large f3, different from 
inclusive DIS, but in agreement with data. 

It is encouraging that simple models reproduce the gross features of the data. 
Given the differences of the models as far as the proton is concerned, it seems that 
hard diffraction probes the wave-function of a virtual photon rather than the structure 
of the proton! 

Hard diffraction in hadron-hadron collisions is much harder to describe and more 
varied, as there can be rapidity gaps between jets, between a jet and a hadron remnant 
Factorization does not seem to hold in this case ||100||, neither is it expected to 



etc.. 



[ p9[| , since, for example, an elastically scattered hadron must traverse the remnant 
of the other hadron, which can cause its break-up. I would like to note, however, a 
recent suggestion [|101|] to describe rapidity gap-like events (between jets) in terms of 
small energy flow in the gap rather than the absence of particles. Although this does 
not correspond exactly to the notion of hard diffractive scattering, such a deflnition 
is more appropriate for a partonic interpretation. 
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4.2 Skewed processes 

Factorization has also been shown for deeply virtual Compton scattering 7*p —>■ 
7P [|102|| and diffractive vector meson production [|103|] (after earlier work in [|104|| ) 
l{Q)p ~* where V can be an onium and Q arbitrary or V can be a longitudinally 
polarized light vector meson, in which case must be large. Note that two-gluon 
exchange is applicable to diffractive vector meson production, but not to diffractive 
DIS, because convolution with the virtual photon wave-function relevant to longitu- 
dinal vector meson production suppresses the asymmetric qq fluctuations, which have 
large transverse size. As a consequence, only the small size qq component contributes 
at leading power. 

Deeply virtual Compton scattering and diffractive vector meson production re- 
quire a generalized parton distribution on the amplitude level, since the proton is 
scattered with non-zero momentum transfer, owing to the difference in invariant mass 
of the initial and final vector particle. These objects, defined as 

j ^e'^'^"- (p'|^(0)7+^(^-)b) (16) 

for quarks, are referred to as skewed (off-diagonal, non- forward, ...) parton distri- 
butions, and describe a parton i (a quark above) extracted from the proton with 
momentum fraction x and returned with momentum fraction x'. For p' = p, the 
skewed parton distribution reduces to the conventional one. The first moment, how- 
ever, is related to a proton form factor. These hybrid properties are also reflected 
in the evolution properties. For x' > the evolution resembles DGLAP evolution. 
For x' < 0, the skewed parton distribution describes emission of a qq pair and the 



evolution resembles ERBL | 105 | evolution of light cone distribution amplitudes. The 
evolution properties and the form of skewed parton densities have been actively stud- 
ied. An interesting observation is that the skewed parton density is determined by 



the conventional one for small x and x' — x ||106 |. 

Is there experimental evidence for skewedness? ZEUS ||107|| reports first evidence 
for deeply virtual Compton scattering, but the data are not yet good enough to allow 
detailed tests. Skewedness effects in diffractive vector meson production are largest if 
the invariant mass difference between the incoming photon and outgoing vector meson 
is large. This suggests to look at T photoproduction ||108|| or vector meson production 
at large ||109|| . Incorporation of skewedness improves the theoretical prediction in 



comparison with data, but other theoretical uncertainties remain large and preclude 
an unambiguous statement. The power behaviour of longitudinal to transverse p 
meson production appears to disagree with the naive estimate ai/ax ~ Q^, but 
the (formally logarithmic) scale-dependence of the gluon distribution may play an 
important role |[109|| . 
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4.3 Exclusive B decays 

There exists a standard framework to discuss exclusive processes at large momentum 
transfer in terms of light cone distribution amplitudes [|105|] . As we are entering the 



era of exclusive B decays, it is only appropriate to consider them as bona fide hard 
reactions. After all they involve momentum transfers ~ ml ~ 25 GeV^ and there 
will be millions of -B's! A general two-body decay amplitude can be written as 

A{B ^ M1M2) = Aie'^'e'^"^' + A2e'^^e'^'^\ (17) 

where Ai^2 denote the magnitudes of the amplitudes, 81^2 strong interaction phases 
and 5wi,2 weak phases. The weak phases are CP- violating and of primary interest. 
Yet to determine them, the strong phases and amplitudes must be known, unless we 
are fortunate enough that only a single term contributes on the right hand side of 
(plTl), or we have enough experimental information to reduce strong interaction input. 

The standard formalism does not immediately apply to B decays, because the 
B meson contains a soft spectator quark. The spectator quark may go to a final 
state meson without participating in a hard scattering. Hence the process cannot 
be described in terms of light cone distribution amplitudes alone. A more general 
factorization for decays into a heavy {D) and a light meson has been proposed in ||110|| , 



based on the idea that the light meson is initially ejected as a compact object from 
the weak decay vertex, although no quantitative conclusions have been drawn, with 
the exception of [|1 1 1|| . Recently, a systematic investigation of the heavy quark limit 



has been undertaken ||112|| . The conclusion is that all soft and coUinear configurations 
can be absorbed either into light cone distribution amplitudes or a form factor for the 
transition B Mi, where Mi is the meson that picks up the light spectator quark. 
In particular, the corrections conventionally termed "non-factorizable" are dominated 
by hard gluon exchange and hence computable. The proposed factorization theorem, 
applicable to heavy-light final states (-Dvr, etc.) and light-light final states (vrvr, ttK, 
etc.), reads 



± 

^(fi^MiMs) = Fb^Mi(O) jdxT\x)^Mix) 



1 







+ / d^dxdyT'\^, X, y) ^B{O^Miy)^Mix) 







with corrections that are suppressed as AQci)/^b- The second term is present only 
for light-light final states and has the form of a standard BL-type term. It accounts 
for hard gluon interactions with the spectator quark ||113 . 



The implications of (0) taken at face value are far-reaching. Since non-pertur- 
bative form factors and light cone distribution amplitudes can either be measured 
or determined in principle with lattice QCD, the strong phases 61^2 and amplitudes 
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Ai^2 are completely predicted. CKM parameters can then be directly extracted from 
measurements of branching fractions and CP asymmetries. 

Some work remains to be done to demonstrate that (p!8D gives accurate predictions 
at the b quark scale. A factorization proof to all orders has yet to be given, which 
may imply that an integration over intrinsic transverse momentum in the B meson 
has to be added to the second term of (pIS]) . Power corrections in AqcB/^b can turn 
out uncomfortably large, if enhanced by small current quark masses. It is also worth 
noting that the light cone properties of B mesons have remained largely unexplored. 
In any event, the new approach improves over naive factorization, which has been 
the most commonly used theoretical tool. Because of its potential for B factories, 
applications need to be carefully examined. 

Conclusion 

QCD is a lively field of incredible variety. It is also often technical. Comparing today's 
QCD overviews with the discussion of big ideas 25 years ago, this variety may even 
appear intimidating. But this transformation in style reflects like no other indicator 
the progress in understanding how QCD works. The challenges provided by strong 
coupling have led to insights into how field theory works unparalleled by any other 
theory. Given the intrinsic beauty and simplicity of QCD, together with its role in 
the future high energy physics programme at the energy frontier, we can be sure of 
further progress in the field. 

I could not have given this talk without the help of many collegues. I benefitted in 
particular from discussions with G. Buchalla, S. Catani, J. Forshaw, A. Hebecker, 
M. Kramer, Z. Kunzst, G. Salam and T. Teubner. 
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Discussion 



Giinter Grindhammer (MPI, Munich): Considering deep inelastic scattering at 
low X, what happens with the difference between the BFKL and DGLAP approaches 
in the case of heavy quark production? In this case one has two hard scales, and 
the quark mass. 



The answer depends on where the heavy quarks are. If the heavy quarks couple to 
the virtual photon ("the top of the ladder"), there is no difference to inclusive DIS 
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at small x. If the heavy quark pair is coupled to the bottom of the ladder, both ends 
are perturbative, as is the case for forward jet production at pt S> Aqcd- In this case 
there is no reason to believe that DGLAP evolution should be relevant. 
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